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Middle Eastern respiratory syndrome, an emerging viral infection with a global case fatality 
rate of 35.5%, caused major outbreaks first in 2012 and 2015, though new cases are continu- 
ously reported around the world. ‘Transmission is believed to mainly occur in healthcare set- 
tings through aerosolized particles. This study uses Quantitative Microbial Risk Assessment 
to develop a generalizable model that can assist with interpreting reported outbreak data or 
predict risk of infection with or without the recommended strategies. The exposure scenario 
includes a single index patient emitting virus-containing aerosols into the air by coughing, 
leading to short- and long-range airborne exposures for other patients in the same room, 
nurses, healthcare workers, and family visitors. Aerosol transport modeling was coupled with 
Monte Carlo simulation to evaluate the risk of MERS illness for the exposed population. 
Results from a typical scenario show the daily mean risk of infection to be the highest for 
the nurses and healthcare workers (8.49 x 10-4 and 7.91 x 10~*, respectively), and the low- 
est for family visitors and patients staying in the same room (3.12 x 104 and 1.29 x 10°, 
respectively). Sensitivity analysis indicates that more than 90% of the uncertainty in the risk 
characterization is due to the viral concentration in saliva. Assessment of risk interventions 
showed that respiratory masks were found to have a greater effect in reducing the risks for all 
the groups evaluated (>90% risk reduction), while increasing the air exchange was effective 
for the other patients in the same room only (up to 58% risk reduction). 
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1.1. Historical Background 


Coronaviruses (CoVs) are a common cause 
of upper respiratory infections in humans. Strains 
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endemic to the human population include 229E, 
HKU, NL63, and OC43 (Chan et al., 2015), which 
circulate continuously among children and adults 
worldwide with infection trends typified by season- 
ality and mild symptoms in healthy individuals. The 
new millennium, however, has marked the emer- 
gence via zoonosis of two highly virulent CoV strains 
novel to the human population. In 2003, a novel CoV 
emerged in the Guangdong Province of China that 
caused a new and deadly outbreak of respiratory 
disease in humans termed as severe acute respira- 
tory syndrome (SARS-CoV). Within months, SARS- 
CoV spread rapidly to 25 countries (in part due to 
the modern, highly globalized nature of air travel), 
with thousands sickened and close to 800 fatal cases 
(Hilgenfeld & Peiris, 2013). Investigation of the ori- 
gin of SARS-CoV led first to the exotic animal mar- 
kets of China and the initial implication of palm civet 
cats and raccoon dogs (which were found to be inter- 
mediate hosts), with further study indicating bats as 
the true natural reservoir of SARS-like CoVs (Han 
et al., 2015). In April 2012, an outbreak of severe res- 
piratory viral illnesses localized in several intensive 
care units occurred in the Middle Eastern country of 
Jordan; both patients and healthcare workers were 
infected. Within several months, cases had also sur- 
faced in several nearby Middle Eastern countries in- 
cluding Saudi Arabia, Qatar, and the United Arab 
Emirates, with rapid movement into 20 additional 
countries in North Africa and Europe (AI-Tawfiq, 
2013; Breban, Riou, & Fontanet, 2013). Early indi- 
cations pointed to a highly virulent infectious agent, 
as a high percentage of patients were dying, partic- 
ularly those with comorbidities. A novel strain of 
CoV was soon isolated and named after the region of 
origin in conjunction with the primary manifestation 
of symptoms—Miuddle Eastern respiratory syndrome 
CoV (MERS-CoV) (de Groot et al., 2013). 


1.2. MERS-Related Health Issues 


MERS-CoV affects the lungs and respiratory 
system with an estimated 35.5% mortality in patients 
globally (World Health Organization, 2018). There 
are currently no human vaccines available to counter 
infection with MERS-CoV, while veterinary vac- 
cines for camels are currently under developments 
(Widagdo, Okba, Stalin Raj, & Haagmans, 2017). 
Therefore, to date, containment of infectious viruses 
via personal hygiene, use of personal protective 
equipment (PPE), isolation of MERS symptomatic 
persons, and quarantine of potentially exposed 
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individuals to prevent contact with others is rec- 
ommended (CDC, 2017a, 2017b). The published 
literature on MERS has consistently estimated a 
reproductive number (Ro; the average number of 
secondary cases generated by a primary case) of <1, 
suggesting that MERS-CoV does not yet pose a pan- 
demic risk (Breban et al., 2013; Nishiura, Miyamatsu, 
Chowell, & Saitoh, 2015; Poletto, Pelat, Lévy-Bruhl, 
Boelle, & Colizza, 2016; World Health Organization, 
2018). In Jeddah, Saudi Arabia, 82 of 168 clinical 
samples stemmed from a single hospital, and phy- 
logenetic analyses of seven confirmed MERS-CoV 
isolates from those cases were found to cluster in a 
single monophyletic clade (Drosten et al., 2015). 

MERS-CoV is primarily transmitted through 
infectious aerosolized particles. Under hospital 
settings, the attack rate has been reported to be 
1.1-10% (Al-Abdallar, 2014; Al-Tawfiq & Perl, 
2015), while 3.6-5% attack rates have been reported 
for the persons in close contact with infected pa- 
tients (Al-Tawfigq & Perl, 2015; Memish, Assiri, & 
Al-Tawfig, 2014). Mean incubation period for the 
virus has been reported to range from 2 to 15 days, 
with a median value of five days (Banik, Khandaker, 
& Rashid, 2015). MERS-CoV infection results in 
fever, cough sore throat, headache, and occasionally 
results in nausea, vomiting, and diarrhea. In more 
severe cases, patients may experience shortness of 
breath, pneumonia, and death (Banik et al., 2015). 
In the South Korean outbreak, the morbidity rate 
was estimated to be 1.08% (Ki, 2015). The patient 
mortality rate has been reported to vary greatly 
depending on the age and underlying conditions, 
such as diabetes, heart disease, and chronic lung 
disease. In the South Korean outbreak, the overall 
mortality rate was reported to be 19.4%. MERS 
infected persons who were already hospitalized for 
other medical conditions had a higher mortality 
rate (33.8%) than the persons without prior medical 
conditions (9.2%). Similarly, patients over 60 years 
of age had a higher mortality rate (38.1%) than 
younger patients (6.4% ) (Ki, 2015). 


1.3. The Large Outbreak in South Korea, 2015 


A cluster of MERS-CoV cases arose in South 
Korea during May 2015. The visitation of a single 
index patient to five different hospitals 1s believed to 
have resulted in 185 downstream nosocomial cases of 
MERS-CoV (Cowling et al., 2015; Park et al., 2015; 
World Health Organization, 2015), although con- 
firmatory phylogenetic analyses have not yet been 
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performed. Unlike the previously documented case 
clusters, the South Korean outbreak was well doc- 
umented with regard to incubation time, transmis- 
sion chains (i.e., 28 first-generation cases, 125 second- 
generation cases, and 32 third-generation cases), and 
contact tracing of infected patients (Ki, 2015). The 
majority of infections were hospital-acquired; only 
one of the 186 patients in the South Korean clus- 
ter was believed to be infected outside of a hospital, 
and two other individuals were infected by modes of 
transmission that are currently unknown (Ki, 2015). 

Despite the fact that MERS has been reported 
to survive a maximum of 24-48 hours on surfaces 
(Van Doremalen & Munster, 2015), it has been 
proposed that based on the South Korean MERS 
outbreak, the virus would not survive long enough 
to be capable of involving spread through indirect 
fomite route (Cho et al., 2016). On contrast, studies 
suggested that the main transmission route of MERS 
was via the airborne route, especially over close 
contact airborne exposure (Xia et al., 2014). Hence, 
isolation of index patient in a negative-pressure 
room and quarantine of potentially exposed persons 
are considered key risk management measures for 
literature that investigated the South Korean MERS 
outbreak (Cowling et al., 2015; Kim et al., 2015; 
Park et al., 2017; Park et al., 2016; World Health 
Organization, 2015). In consequence, isolation and 
quarantine would be measures that would drastically 
lower the risk of MERS infection once patients 
are identified. From previous outbreaks, the index 
patient stays unidentified as a MERS carrier for up 
to two days (Cho et al., 2016). Additionally, the time 
for identifying MERS from a diagnostic laboratory 
in a patient takes up to three days (Cowling et al., 
2015), so probable exposure durations around two to 
three days are relevant scenarios to model. 


1.4. Study Objectives 


The objective of this study was to use the Quan- 
titative Microbial Risk Assessment (QMRA) ap- 
proach to develop a generalizable model for quanti- 
fying the risk of infection associated with in-hospital 
exposures to MERS through infectious aerosols. 
The parameter values were selected from multiple 
sources including the latest reported large outbreak 
that occurred in South Korea, and data from other 
sources. Risk of infection is estimated for four types 
of at-risk populations: nurses and healthcare work- 
ers visiting the index patient (before the patient 
was identified as carrying MERS) and other patients 


sharing the same room; family visitors coming to visit 
the index patient; and the other patients sharing the 
same room (Cho et al., 2016). Risk estimation is con- 
ducted by using the Monte Carlo simulation method 
to incorporate uncertainty and variability in the risk 
characterization. Sensitivity of the model parameters 
is assessed to determine where additional data or 
knowledge could potentially reduce uncertainty and 
increase our understanding of these risks. Finally, the 
effectiveness of mask and increased ventilation risk 
management measures is evaluated. Rather than a 
retrospective case analysis, the study is intended to 
contribute a framework for analyzing current and 
future MERS risk in similar settings. 


2 MATERIALS AND METHODS 


2.1. Exposure Scenario and Assessment 


The basis of the exposure scenario involves a 
symptomatic patient infected with MERS-CoV who 
has been admitted to a hospital without implemen- 
tation of isolation or quarantine procedures. It was 
assumed that all exposed people were susceptible to 
infection and all infections led to illness (or death). 
A typical size of 230 m° hospital room was set for the 
model, which is four times the single patient room 
size noted in Yin, Gupta, Zhang, Liu, & Chen (2011) 
and is based on the fact that over 50% of the hospital 
rooms in South Korea have four or more beds. The 
symptomatic patient was considered the only source 
of infection within the room (see Fig. 1). 

MERS-CoV is thought to be transmitted primar- 
ily via aerosols in a manner similar to endemic hu- 
man respiratory CoV strains such as 229E and OC43. 
For the present assessment scenario, only the risk of 
infection from aerosolized particles and droplets ex- 
pelled by coughing was considered. The influence of 
nebulizer treatments that can be done on the index 
patient was considered negligible and not included in 
the model. Although the contribution of this treat- 
ment was suggested by Park et al. (2016), studies 
have also demonstrated that nebulizers do not specif- 
ically impact transmission (Seto, 2015; Thompson 
et al., 2013). Fomites may also serve as a poten- 
tial reservoir for MERS-CoV due to the settling of 
aerosols after release from infected persons. How- 
ever, some studies stated that fomite-based exposure 
pathways were not significant compared to airborne 
routes, and so it was not considered in this study 
(Xiao, Li, Sung, Wei, & Yang, 2018). 
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Strain: SARS 
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Risk —— Outbreak scenario: 2 days, 3 ACH, highest MERS conc. 
is | 
characterization | Risk of infection 
(Pings = 1 — ef), Pap = 1 — (1-Pingh") 
{for nurses, HCVWs, family, potients) 
Wearing of N95 mask —= —— Increased ACH (9, 12 ACH) 
Risk J Risk after intervention 
Management | (Pingia = 1— el", Py = 1 —(1-Ping)”) 


{for nurses, HCW, family, potients) 


Fig. 1. Exposure scenario and QMRA outline steps. QMRA = Quantitative Microbial Risk Assessment; MERS = Middle Eastern respi- 
ratory syndrome virus; HCW = healthcare worker. 
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Two forms of modeling were included in this 
MERS assessment: (1) modeling aerosol concen- 
trations to identify at-risk populations in hospital 
settings; and (2) estimating exposure dose and char- 
acterizing risk. The risk of infection for several ex- 
posure populations was considered as follows: (1) 
other patients in the same room of index patient; (2) 
nurses; (3) other healthcare workers (e.g., doctors) 
visiting the index patient and others in the room; and 
(4) family members coming to visit the index patient. 

Viruses released via coughing and transport in 
the hospital room were modeled using a mass bal- 
ance approach to approximate a steady-state concen- 
tration of viruses contained in aerosol droplets. The 
droplets are being removed from the system either 
due to settling to the floor or ventilation-based air 
exchange. The risk of infection for each of the four 
populations was assessed based on exposures occur- 
ring over 1, 8, 20, and 41 days. These time periods 
were based on reported durations from the symptom 
onset to discharge from the hospital during the Ko- 
rean outbreak—a median of 20 days, minimum of 
8 days, and maximum of 41 days (Ki, 2015)— 
and from estimated durations for other patient 
exposure—up to 44 hours (Cho et al., 2016). 


2.2. Aerosol Transport Modeling 


Aerosol transport modeling was undertaken to 
assess virus inputs from coughing and removal via 
settling onto surfaces and the air exchange pro- 
cesses (i.e., heating, ventilation, and air conditioning 
[HVAC] systems). The model room system was as- 
sumed to have reached steady state, meaning that 
there is no accumulation or loss from the system over 
time, and that the input flow rates must equal the 
removal flow rates. This input—output relationship is 
shown in Equation (1): 


Nin, coughing) — Nout, setting) i Nout, ventilation) 
+N (inhalation) » (1) 


where N is the number of droplets containing viruses. 
In Equation (1), Minhalation) or the number of viruses 
removed through inhalation by infected or unin- 
fected persons (patients in the same room, health 
care workers, and visitors) was assumed to be non- 
significant as compared to the other two terms, 
Nout, settling) and Noout, ventilation) > and thus was ne- 
glected as a sink. Expiratory events (i.e., cough- 
ing) produces a broad distribution of aerosol par- 
ticles, however, this analysis was only concerned 


with aerosols that were likely to be inspirable and 
respirable. Aerosol production values were taken 
from Stilianakis and Drossinos (2010) and the refer- 
ences therein. Particles with a diameter of <10 um 
were considered as respirable aerosols. Respirable 
aerosols are expected to be easily transported, due 
to their small diameter, and thus represent a po- 
tential exposure pathway for people that are farther 
away from the source (e.g., more than 1-2 m from 
the source). Thus, respirable particles were the only 
evaluated exposure pathway for patients sharing a 
room with an infected symptomatic patient. Aerosols 
with a diameter of 10-100 um were considered as 
inspirable aerosols as these large particles are not 
expected to be transported long distances and are 
only relevant for persons in close contact. Nurses, 
healthcare workers, and visitors were assumed to be 
exposed to both respirable and inspirable aerosols. 
Viral release into the room was calculated using 
Equation (2): 
md? 


V; = z x 107", (2) 


where V; (mL) is the volume for each droplet size d; 
that are released into the room as inspirable or res- 
pirable droplets during each coughing event. Each 
cough produced N; number of droplets of size d;, 
where each droplet is assumed to be spherical, and 
the droplet volume is calculated as ; x a x d?, where 
d; is the diameter (um). The droplets were assumed 
to be produced from a patient lying supine, such that 
the droplet cloud was produced at a 1 m height. 

Following Stilianakis and Drossinos (2010), 
pathogen generation (e.g., coughing) and removal 
(e.g., settling, ventilation) were assumed to be a 
continuous process. Exhalation by the infected pa- 
tient was not considered a source of virus-containing 
droplets. 

After the particles were produced during a 
coughing event, droplet evaporation, droplet set- 
tling, and droplet removal via the ventilation were 
considered. Postevaporation particle transport was 
evaluated, accounting for two removal mechanisms: 
droplet settling and ventilation-based droplet re- 
moval. Stoke’s law was used to calculate droplet 
terminal settling velocity vj(terminal) (m/hr) (Equa- 
tion (3)), which was assumed to be impacted only 
by particle diameter d; (Nicas, Nazaroff, & Hubbard, 


2005). 
0.166 
di | ©) 


Vi(terminal) = 0.108 x dx [14 


Terminal settling velocities were calculated for each 
of the representative particle sizes, d;. A critical set- 
tling velocity, v¢criticaty (m/hr), was calculated as the 
required settling velocity to fall from the height of 
the patient bed /ycoush) (m) during the air residence 
time t (hour) (Equation (4)). Air residence time, T, 
is the average amount of time that a “parcel” of air 
is in the room, which depends on the volume of the 
TOOM U(yoom) and the ventilation rate (ventilation): 


es U(room) (4) 
(ventilation) 
A(cough 

U(critical) = —_ ; (5) 


Ventilation flow rate ventilation) WaS quantified by 
the number of air exchanges per hour (ACH) of the 
room volume, which was defined as shown in Equa- 
tion (6): 


(ventilation) = Vroom) x ACH, (6) 


where Viroom) is the volume of the hospital room (m*) 
and (exchange) 18 the air flow rate (m°/hr) determined 
by the number of ACH vvexchange). AS Stated previ- 
ously, height of the patient bed, h(cough), Was 1 m. For 
particles that had a terminal settling velocity greater 
than the critical velocity (vj (terminal) > U(critical)), 1t Was 
assumed that settling was a viable removal mecha- 
nism. It was further assumed that droplets that hit 
the floor were permanently removed from the sys- 
tem with no resuspension. This acknowledges that all 
the settleable aerosol droplets settled to the floor in a 
time interval less than t. However, due to the contin- 
uous generation, there were some fractions of the set- 
tleable droplets that were not yet settled. At a given 
time, for the droplets with terminal velocity greater 
than the critical velocity (vj (terminal) > U(critical)), 1 Was 
assumed that the aerosol concentration of settleable 
droplets was proportional to the ratio of settling ve- 
locities, as shown in Equation (7). For the droplets 
that had terminal settling velocities less than the crit- 
ical settling velocity (vj(terminal) < U(critical)), It Was as- 
sumed that there was no droplet removal via settling. 


Ni(room, settleable) — Nicin, cough) | Pee . (7) 

Uj (terminal) 
For these later particles, it was assumed that air cur- 
rents in the room dictated their transport. However, 
this transport and homogeneous mixing did not in- 
clude settling onto another surface resulting in re- 
moval (i.e., striking a piece of furniture, or a wall) 
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and was considered entirely an elastic collision. Air- 
borne particles were assumed to be homogeneously 
distributed within the volume of the room. Hence, 
the number of droplets containing viruses removed 
through settling for each droplet 7 1s: 


U(critical) | . (8) 


Vi (terminal) 


Nicout, settling) = Nicin, cough) 1 = 


Air exchanges via ventilation was also consid- 
ered a removal mechanism, in which air, including 
the homogeneously mixed virus-containing aerosol 
droplets, was removed from the hospital room and 
replaced with new air. It was assumed that the re- 
placement air contained no viruses. During each 
air replacement, all the remaining droplets were as- 
sumed to be removed by the ventilation, which 1m- 
plies the relationship in Equation (9). 


Ni(out, ventilation) = Nicin, cough) —_ Ni(out,settling)« (9) 


It was assumed that Ny(out,ventilation), NUMber of 
droplets remaining after settling, are suspended in 
the room until they are removed by ventilation. 
Hence, the concentration of saliva in the air pro- 
duced by a single cough per unit volume of room air 
is calculated as 


ee 


10 
Vroom) ( 


C(saliva in air) = 
where C,saliva in air) 1S the concentration of saliva in 
the room air produced by a single cough per hour 
(mL/m*), V; is the volume of each droplet calculated 
using Equation (2), and Vizoom) is the room volume. 

We further assumed a standard air exchange rate 
of six times per hour (Zumla & Hui, 2014). The half- 
life of CoVs in the air is 67.33 hours (Ijaz, Brunner, 
Sattar, Nair, & Johnson-Lussenburg, 1985), but since 
we assumed that the air in the room was exchanged 
six times per hour, decay was not considered. 


2.3. Aerosol Concentrations in the Air 


To model the amount of virus released into the 
air, several studies were compared that specified 
the number and size of droplets expelled during 
coughing (Duguid, 1946; Loudon & Brown, 1967; 
Nicas et al., 2005; Papineni & Rosenthal, 1997) for 
selection of the data set that best fits the condition 
of patients exposed to MERS-CoV. The number 
of cough events per hour was modeled based on 
Loudon and Brown (1967), using the estimates for 
the number of cough events in nonsmokers with 
pneumonia. Based on Nicas and Jones (2009), we 
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assumed that 0.044 mL of saliva was emitted per 
cough, which represents the most conservative es- 
timate compared to other published volumes in the 
literature (Duguid, 1946; Loudon & Brown, 1967; 
Papineni & Rosenthal, 1997). Saliva volume was 
assumed to have a uniform distribution with a +10% 
of the reported value. Of the expired fluid, 0.00015 % 
was considered respirable and 0.54% was considered 
inspirable. In other words, about 99.45% of the 
volume expired during each cough was considered to 
be nonrespirable and noninspirable, and therefore 
was not included in this analysis. Respirable droplets 
were modeled as aerosols with mean postevapora- 
tion diameters of 4 um and 8 um (for small and large 
respirable droplets), which Stilianakis and Drossinos 
(2010) estimated were produced at a rate of 160 
and 7.5 droplets per coughing event, respectively. 
Similarly, based on Stilianakis and Drossinos (2010), 
representative inspirable droplets corresponded to 
aerosols with mean postevaporation diameters of 
7.3 um and 74 um diameter droplets (corresponding 
to inhalable aerosols), which were produced at 41.47 
and 138.48 droplets per cough, respectively. Other 
than this initial evaporation, it 1s assumed that the 
aerosol droplets did not change in size, including that 
neither further evaporation nor particle aggregation 
occurred. Uncertainty in the droplet production 
numbers was investigated by holding the number 
of particles constant, and using bootstrap iterations 
to compare the uncertainty in the relative number 
of particles for each of the four respective repre- 
sentative particle sizes. The results of the bootstrap 
uncertainty analysis were used to model particle 
production as a stochastic input. 


2.4. MERS-CoV Concentration in Saliva 


Multiple papers have quantified levels of MERS- 
CoV in sputum, nasopharyngeal secretions, and 
saliva samples using the quantitative polymerase 
chain reaction (qPCR) methodology (Corman et al., 
2015; Min et al., 2016; Muth et al., 2015). MERS- 
CoV titer data specified in these studies are in to- 
tal viral units (noninfectious + infectious) of RNA 
genomic copies per milliliter (GC/mL) as the val- 
ues were generated using real-time qPCR. Since the 
dose-response model unit was in plaque-forming 
unit (PFU), according to the used best-fit dose— 
response for SARS-CoV taken from the QMRA 
Wiki website (Huang, 2013), a conversion factor of 
1,239:1 (1,239 GC equivalent units to one infectious 


PFU) reported by Houng et al. (2004) and based on 
a SARS-CoV qPCR assay was employed to calcu- 
late infectious PFU values for the MERS-CoV ex- 
posure modeling. Recovered MERS-CoV concentra- 
tion data were fitted to a lognormal distribution. 


2.5. Exposed Population Behavior 


Exposure scenarios for the nurses and healthcare 
workers were modeled based on the frequency and 
duration of their patient visits. For healthcare per- 
sonnel, due to the wide range of reported durations 
per visit by Cohen, Hyman, Rosenberg, and Larson 
(2012), a triangular distribution was specified with a 
median of two minutes and a range of 1-72 minutes) 
(Table I). Similarly, a triangular distribution with a 
median value of two minutes per visit and a range 
of 1-120 minutes was assumed for the nurses as in- 
puts in the exposure model (Cohen et al., 2012). For 
both the healthcare workers and nurses, the number 
of patient visits and number of different patients vis- 
ited were also taken from Cohen et al. (2012) and 
are tabulated in Table I with all model inputs and 
distributions. Nurses and healthcare workers were 
assumed to be exposed to inspirable and respirable 
particles while visiting the index case, and to the res- 
pirable particles while visiting other patients in the 
same room. Other patients in the room were assumed 
to be exposed to respirable particles only 24 hours 
a day (Ki, 2015). For the family visitors, a median 
visit duration of 14 minutes was used (Cohen et al., 
2012). Furthermore, based on Cohen et al. (2012), 
frequency of visitors was assumed to range from 0 to 
6.4 visits per hour with a median value of 1.3. Daily 
exposure doses for nurses, healthcare workers, the 
other patients, and family visitors were calculated by 
ageregating the exposure doses over the entire day 
consisting of multiple visits. 


2.6. Estimated Exposure Dose 


The daily exposure dose for the nurses and 
healthcare workers was calculated by considering 
that once entering the room, they would expose 
themselves both through respirable and inhalable 
aerosols during their visit to the MERS index pa- 
tient, and through only respirable aerosols when 
visiting the other patients in the room. Hence, daily 
exposure dose for nurse and healthcare worker 
consisted of the sum of each of these two exposure 
routes: 
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Table I. Parameters Used in the Model 

Parameters Unit Description Input Values (a; b)* Distribution Sources 

Vsaliva/cough) mL Volume of saliva 0.044 (0.0396; 0.0484) Uniform Nicas and Jones 
expelled/cough (+10%) (2009) 

R(GC:PFU) _ Genomic copies-to-PFU b2o971 Point value Houng (2004) 
conversion factor 

C(MERS in saliva) PFU/mL _ Virus conc. saliva = Conc. 41,734 (7; 201,945) Lognormal Corman (2015), 
[#GC/mL] x R(GC:PFU) Min (2016), 

Muth (2015) 

N(cough/day) day! Number of coughs/day = 625 (0.125; 39.25) Triangular Loudon (1967) 
N(cough/hr) x 24 

d; um Droplet diameter (4 um 4; 8; 7.3; 74 Point value Stilianakis and 
and 8 um for <10 um Drossinos (2010) 
respirable droplets, 

7.3 um and 74 um for 10- 
100 um inspirable) 

N; # Number of 160; 7.5; 41.47; 138.48 Point value Stilianakis and 
droplets/diameter d; Drossinos (2010) 
emitted/cough 

V; mL Volume of each Calculated Point value Stilianakis and 
droplet/diameter d; = Drossinos (2010) 
(x d;*)/6 x 10-7 

Vou) m Hospital room size 230 Point value —- Yin (2011) 

U(critical) m/hr Required droplet settling Calculated Point value Nicas (2005) 
velocity to fall on ground 
= 0,108 x4? 

(1 + 0.166/d;) 
C(saliva in air) #/m?> Conc. droplets in the Calculated Normal Stilianakis and 
air/cough Drossinos (2010) 
n 
= (2 Ni Vi)/Viroom) 
= 

N (room entries/hr) he Visit frequency of nurse 2.5 (0512.0) Triangular Cohen (2012) 

N room entries/hr) he Visit frequency of 1.6 (073.12) Triangular Cohen (2012) 
healthcare workers 

N (room entries/hr) hr! Visit frequency of a family 1.3 (0; 6.4) Triangular Cohen (2012) 
member 

N(patients visited /entry) # Number of different 4.5 (0.5; 18) Triangular Cohen (2012) 
patients visited by a 
nurse 

N(patients visited /entry) # Number of patients visited 2.8 (0.5; 7) Triangular Cohen (2012) 
by a healthcare worker 

t(spent/entry) min Time spend/Visit of a nurse 2 (1; 120) Triangular Cohen (2012) 

(spent /entry) min Time spend/visit of a 3 (1; 72) Triangular Cohen (2012) 
healthcare worker 

t(spent/visit) min Time spend/visit of a family 14 (1; 124) Triangular Cohen (2012) 
member 

texposed/d hr/d Contact time of other 24 Point value Assumed 
patient in the same 
room/d 

Viinhaled/d) m?/hr Respiration rate of an 0.5 Point value EPA (2011) 
exposed person 

k PFU"! Parameter of the 0.00246 (0.00135; 0.00459) Normal Huang (2013) 
exponential 
dose-response 

ACH hie Air exchange rate (for the 6 Base case Zumla and Hui (2014) 
base scenario) 

F(aroplets out mask) % % droplets out mask (from 0.032 (0.010; 0.100) Uniform Borkow (2010), Wen 


log reduction) 


(2013) 


“a = Min value for triangular and lognormal distribution and 5th percentile value for normal distribution, respectively; b = max value for 
triangular and lognormal distribution and 95th percentile value for normal distribution, respectively. 
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Dexpo/d,n—hew) = C(MERS in saliva) X C(saliva in air) 


x Viinhaled/d) 


x Ncough/hr) x 
(ventilation) 


N (room entries/hr) 
x — ss X U(gpent/entry) X L(work/d)> 


(11) 


where Deexpo) 1S the daily MERS virus inhaled by 
exposed personnel while being one time near index 
patient and another time near patients sharing the 
room (PFU/day), Camers in saliva) 18 the concentration 
of MERS in saliva (PFU/mL), C(sativa in air) 18 the 
concentration of droplets in the air after one cough, 
Ncoughs/hr) 18 the number of coughs per hour (#/hr), 
(ventilation) 18 the ventilation air flow rate of the 
room derived from the ACH (#/hr), Vinhalea/d) 1S 
the air intake rate of the exposed person (m*/hr), 
Noom entries/hr) 18 the number of entries nurse or 
healthcare worker makes per hour to visit either 
the index patient or the other patients (#/hr), 
N(patients visited/entry) 1S the number of patients visited 
by nurses or healthcare workers per room entry visit 
(for the index patient or other patients) (#/visit), 
t(spent/entry) 1S the amount of time spent during each 
visit (hr/visit), and f(work/a) 18 the number of daily 
working hours for nurses and healthcare personnel 
(assumed 8 hr/day). 

For the other patients in the same room, daily 
exposure dose was calculated as follows: 


N(patients visited/entry) 


Daexpo,op = C(MERS in saliva) * C(saliva in air) X N cough /hr) 


1 
ee V(inhaled/d) x lexposed/d> 
(ventilation) 


(12) 


where the daily exposure duration fexposeasd WaS as- 
sumed to be continuous (i.e., 24 hr/d). 

For the family visitors, daily exposure dose was 
calculated based on their number of visits per day of 
the index patient Namily visits/d: 


Dexpo,fim = C(MERS in saliva) * C(saliva in air) X Ncough/hr) 


1 
x ———_ X V(inhaled/d) x N family visits/d) 
(ventilation) 


a (spent /visit) . (13) 


A systematic literature review was conducted to 
determine the best estimates for each input param- 
eter in the exposure model. A Monte Carlo simula- 


tion was conducted using the Crystal Ball® program 
(Version 11.1.4512.0, Oracle, Redwood Shores, CA, 
USA) to incorporate variability and uncertainty in 
the input parameters and to propagate it to the out- 
put parameters (i.e., exposure doses per subpopula- 
tion, risks of infection). Risks of infection for each 
scenario were calculated using a published dose— 
response model as described in Section 3. A differ- 
ential sensitivity analysis of model variance was per- 
formed to determine which input variables have the 
ereatest effect on the risk estimates. To reduce the 
risk of MERS infection, two types of risk mitigation 
strategies were evaluated using the final risk models: 
increasing air exchange rate and using a mask as PPE. 


2.7. Dose-Response Model 


A primary knowledge gap in the study is the 
absence of a dose-response model for MERS-CoV. 
Therefore, the SARS dose-response model (Huang, 
2013) was employed as a surrogate. MERS has sev- 
eral similarities to SARS: both have an animal ori- 
gin and appeared around 2002 in approximately 
the same regions—Asia and Middle East (Sutton 
& Subbarao, 2015), both are respiratory CoVs with 
the same transmission route, both have a compara- 
ble protein structure for binding to host cells (Lu, 
Wang, & Gao, 2015), and both have reported simi- 
lar tropism within cells (Zhou, Chu, Chan, & Yuen, 
2015). Hence, despite probable differences in attack 
rates and mortality rates between the two viruses 
(Chan et al., 2015), it was assumed in this study that 
the SARS dose-response model is the best avail- 
able model for MERS. Several dose-response stud- 
ies for SARS were evaluated to determine a rec- 
ommended dose-response model (De Albuquerque 
et al., 2006; DeDiego et al., 2008; Mitchell & Weir, 
n.d.; Watanabe, Bartrand, Weir, Omura, & Haas, 
2010). Recommended SARS dose-response model 
follows the exponential dose-response relationship 
(Equation (14)) for exposure dose expressed in PFU 
and the probability of a response based on an end 
point of death in mice (De Albuquerque et al., 
2006; DeDiego et al., 2008). For translating this an- 
imal dose-response relationship to a human dose— 
response relationship, a generally accepted assump- 
tion that a death end point for an animal model may 
be used for examining the human risk of infection 
was applied (Haas, Rose, & Gerba, 2014). 

The general equation for the exponential model 
is: 


Pont = 1-e*, (14) 
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where Pry is the risk (probability) of infection, k 
is the optimized dose-response function parameter 
(PFU~'), and d is the dose (PFU). In the Monte 
Carlo analysis, the k value in the dose-response 
model was modeled with a normal distribution based 
on the 5th, 50th (median), and 95th percentile val- 
ues reported by Huang (2013) and reported in 
Table I. 

The cumulative risk of the morbidity across mul- 
tiple exposure days was modeled by Equation (14) 
(Haas, Rose, & Gerba, 2014): 


Py = 1—(1— Pui)”, (15) 


where Py is the probability of morbidity and n is the 
number of days of exposures with Pip¢, probability 
of infection from a daily exposure. The risk associ- 
ated with each population was assessed for 8, 20, and 
41 days of exposure, which represents the minimum, 
median, and maximum hospitalization periods for an 
MERS infected patients (Ki, 2015). 


2.8. Risk Management Evaluation—Air Change 
Per Hour and Wearing of Mask 


To reduce the amount of airborne respirable par- 
ticles, Zumla and Hui (2014) recommend increasing 
the air changes per hour (ACH) from 6 to 12 in hos- 
pital facilities or rooms with high risk of airborne 
disease. Thus, in addition to the worst-case scenario 
considering 0 ACH and the Korean outbreak sce- 
nario using 3 ACH (Cho et al., 2016), standard 6 
ACH (Zumla & Hui, 2014), along with increased 9 
and 12 ACH were evaluated for their efficacy in min- 
imizing the infection risk. 

We evaluated the use of respiratory masks (N95) 
as a means of personal protection. Laboratory stud- 
ies showed a large decrease (up to >4 log reduction) 
in virus exposure when wearing masks (Borkow, 
Zhou, Page, & Gabbay, 2010). However, the de- 
crease did not take into account imperfect mask fit 
or lack of compliance in wearing the masks. Due 
to these factors, MERS-CoV reduction due to wear- 
ing N95 respirators was assumed to have a uniform 
distribution spanning 1-2 log reductions in MERS- 
CoV concentration (Balazy, Toivola, Adhikari, et al., 
2006; Batazy, Toivola, Reponen, et al., 2006; Gupta, 
2011; Rengasamy, Zhuang, & Berryann, 2004; Wen 
et al., 2013). 
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3. RESULTS 


3.1. Risk of Infection 


Based on the results obtained from the Monte 
Carlo simulation, the risk of infection to each 
exposed group was characterized. Fig. 2 shows the 
boxplot of daily risk of infection to each group for the 
base scenario, meaning without any preventive in- 
terventions, the standard rate for hospitals of 6 ACH 
was considered (Zumla & Hui, 2014). The median 
(mean) daily risk for the nurses coming to visit the 
index patient and other patients in the same rooms; 
the healthcare workers (e.g., doctors); the family 
members coming to visit the index patient; and the 
other patients sharing the room were found to be 
1.33 x 10-® (8.49 x 1074), 1.18 x 107° (7.91 x 1074), 
6.36 x 10~? (3.12 x 10~*), and 2.73 x 10~? (1.29 x 
10-*) respectively. The estimated highest daily risk 
of infection for the healthcare workers and nurses 
suggested the frequency of airborne close-range ex- 
posure route plays a bigger role in the transmission of 
MERS, compared to the long-range airborne route 
to which other patients are exposed, confirming what 
was suggested by Xiao et al.’s (2018) work. Statistical 
t-tests showed that the daily risk of infection for 
healthcare workers was significantly higher than 
the one for the other patients or the family visitors 
(p-value = 0.0014 and 0.0240, respectively, at a = 
0.05). When comparing nurses and other healthcare 
worker, the result is not significant (p-value = 
0.8475), so they have similar risks. Other patients in 
the same room had a statistically significant lower 
risk of infection compared to nurses (p-value = 
0.0017), but had nonsignificant statistical differences 
in risk with family visitors (p-value = 0.0547). 

Fig. 3 shows the aggregated risk of infection 
for the exposed populations during multiple daily 
exposures to the MERS-infected patient for typical 
hospital durations. As expected, the results show 
increased risk of infection to all the exposed popula- 
tions over time. The rate of increase was highest for 
the healthcare workers and nurses, in comparison 
to the family visitors, which itself had higher rate of 
increase compared to the other patients sharing the 
room. Similar to the daily risk, aggregated risk of 
infection was the highest for the healthcare workers 
and nurses, followed by family visitors and other 
patients. By day 41, the average risk of infection to 
the nurses was 1.01, 1.2, and 2.4 times the risk for 
the healthcare workers, family members, and other 
patients, respectively. 
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Fig. 2. Boxplots of daily risk of MERS 
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Fig. 3. Aggregated mean risk of infection rate with six air changes over days exposed. Vertical bars denote standard error. HCW = health- 


care worker. 


3.2. Parameter Sensitivity 


Fig. 4 shows the parameter sensitivity of the 
model for daily risk of infection. A parameter with 
a greater rank correlation coefficient indicates that 
the input parameter distribution was more correlated 
with the output risk of infection for the population 


specified. Input parameter sensitivity can either be 
due to the uncertainty in estimating the value of a 
parameter, the known naturally occurring variance 
of this parameter, or because disparate data from 
different sources were used to estimate the range 
of a parameter. Here, while many parameters were 
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Sensitivity analysis - Input parameter contribution to risk output variance 
(for 1 day, at 6 ACH) 
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Fig. 4. Parameter sensitivity of the daily risk. ACH = air change per hour; HCW = healthcare worker. 


modeled stochastically, only the most sensitive pa- 
rameters are shown in Fig. 4. 

For all exposed groups, the concentration of 
MERS virus in the saliva was the most sensitive pa- 
rameter, which accounted for over 90% of the uncer- 
tainty in the daily risks. Because viral load of MERS 
in saliva is believed to naturally vary among peo- 
ple, sensitivity of this parameter is understood as be- 
ing due to natural variance. For the family members, 
nurses, and the healthcare workers, the duration of 
each visit was the second most sensitive parameter, 
with contribution to variance ranging from 1.7% to 
2.6%. As the other patients were assumed to be con- 
tinuously exposed, this parameter was obviously not 
important for them. 

Other parameters contributing to the risk vari- 
ability were the rate of coughing per hour from the 
infected index patient (about 1.5% contribution), the 
visit frequency of healthcare workers, nurses or fam- 
ily (about 1.2% contribution), and the dose-response 
parameter k (about 0.5% contribution). 


3.3. Risk Management Evaluation 


The results of the risk management evaluation 
showed that increasing the air ventilation rate from 
6 to 9 or 12 ACH was an effective risk mitigation 
measure for the other patients sharing the room with 
the index patient, but not for the other persons in 


close contact (Fig. 5). For the other patients, mean 
daily risk of infection could be reduced by about 
30% or 58% through increasing the air ventilation 
from 6 to 9 or 12 ACH, respectively. For the nurses, 
healthcare workers, and family visitors, only up to 
about 2% reduction in mean daily risk could be 
achieved by increasing the ACH from 6 to 12. 

Using a mask was found to be the most effective 
intervention measure in minimizing the risk of infec- 
tion. By using the mask, about 89-97% of the mean 
daily risk could be reduced for other patients, nurses, 
healthcare workers, and family visitors. Higher risk 
reduction suggests that all the exposed groups should 
use mask as a PPE to minimize the associated risk of 
infection. 


4. DISCUSSION 


With the recent emergence of MERS-CoV on 
the global scene, much remains unknown about the 
way the virus behaves. Animal models are being 
evaluated for their suitability for dose-response 
models, and nonhuman-primates appear promising, 
exhibiting similar symptoms to human (Sutton 
& Subbarao, 2015). However, no dose-response 
models have been completed for MERS. Hence, the 
model developed applied a dose-response relation- 
ship based on SARS-CoV pathogen, believed to be 
the best surrogate to use for MERS-CoV because 
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Comparison of risk reduction per intervention 
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Fig. 5. Effect of risk management strategies on average daily risk of infection. ACH = air change per hour; HCW = healthcare worker. 


of their identical animal origin, same regions of 
main occurrence in Asia and Middle East (Sutton & 
Subbarao, 2015), their consistent respiratory trans- 
mission route, their comparable protein structure for 
binding to cells (Lu et al., 2015), and their similar 
tropism within cells (Zhou et al., 2015). For these 
reasons, the proposed best-fit SARS-CoV dose-— 
response from QMRA Wiki website was applied 
(Huang, 2013). However, it must be mentioned that 
some other researchers applied different SARS-CoV 
dose-response, such as the Xiao et al. (2018) team 
did. In addition, since depending on the context, 
SARS-CoV and MERS-CoV may show different or 
higher infectivity (Chan et al., 2015), the question 
could arise to maybe use a developed dose-response 
from other RNA viruses, such as influenza (Chabre- 
lie, Mitchell, Rose, Charbonneau, & Ishida, 2018) 
or any other. Yet, estimated risk of infection 1s 
believed to slightly vary if using other SARS-CoV 
dose-responses, because of their relatively similar 
k parameter value: Xiao et al. (2018) referring to 
a k value of 0.0032, while DeDiego et al. (2008) 
estimated a median k value of 0.00297 and De Albu- 
querque et al. (2006) estimated a median k value of 
0.00214, all to be compared with the median k value 
of 0.00246 derived from the best-fit SARS-CoV 
dose-response provided by QMRA (Huang, 2013). 
Assumptions about the frequency and length of 
visits to patients in the hospital were based on the 
best available data. The visitor data (i.e., frequency 


and duration) were obtained from a hospital in New 
York City (Cohen et al., 2012), which may not accu- 
rately represent the visitor behavior in all countries. 
For instance, in South Korea the healthcare system 
can give responsibility to family members of the pa- 
tients to provide some care for the patients, leading 
hospital rooms to be often crowded with patients and 
their families or privately hired healthcare aids tak- 
ing care of the patients (Ki, 2015). 

Regarding the RNA copies:PFU ratio parame- 
ter, no ratio was reported from literature for MERS- 
CoV, but several were for SARS-CoV, ranging from 
1:1 (Xiao et al., 2018) to 300:1 (Sampath et al., 2005), 
and up to 1,200-1,600:1 (Houng et al., 2004). In con- 
sequence, this parameter appears to possibly change 
the risk estimation. Running the model with a 1:1 
ratio gave daily risk of infection at 6 ACH about 2 
log higher than the ones reported with a 1,239:1 ra- 
tio (developed base model), with new mean daily risk 
of 3.07 x 10-7, 2.96 x 10~7, 2.26 x 10-7, and 1.39 x 
10-* for nurses, healthcare workers, family visitors, 
and other patients, respectively (data not shown). 
Applying a uniform distribution from 1:1 to 1,600:1 
for the copies:PFU ratio lead this parameter to con- 
tribute for risk variance to about 3%. Because the 
RNA copies:PFU proposed by Sampath et al. (2005) 
was derived from SARS-CoV isolated from multiple 
different animals, and because Xiao et al. (2018) ap- 
plied a 1:1 ratio simply because they recognized no 
ratios were reported for MERS-CoV, this study used 
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the median ratio of 1,239:1 proposed by Houng et al. 
(2004). 

Additionally, only the airborne exposure trans- 
mission route through inhalation was considered in 
the model. Despite the fact that literature suggests 
aerosol-generating procedures such as bronchoscopy 
or intubation might amplify viral transmission (Tran, 
Cimon, Severn, Pessoa-Silva, & Conly, 2012), this 
study did not consider this aspect, as a simplifying as- 
sumption to develop a simple and easily adaptable 
model. Since an index patient in the MERS was re- 
ported having undergone a bronchoscopic examina- 
tion and nebulizer therapy (Park et al., 2016), risk 
of infection might be underestimated. However, it 
must be clarified that according to Seto’s (2015) and 
Thompson et al.’s (2013) works, nebulizer usage have 
been proven to not significantly increase the viral 
transmission through air. Other exposure routes of 
infection, such as the fecal-oral and fomite transfer 
routes that have been postulated for other CoVs, 
were not included in this analysis. This decision was 
motivated by the fact that literature suggests that 
MERS transmission mainly occurs through the short- 
range airborne exposure route (Xiao et al., 2018). 

To assess the validity of the model developed, 
a scenario mimicking the South Korean conditions 
was tested. The outbreak selected was one of the 
two major outbreaks that occurred in South Korea, 
specifically the one reported in the Samsung Medical 
Center hospital outbreak in Seoul, for which attack 
rates per population exposed were given from litera- 
ture (Cho et al., 2016). No viral loads were reported 
from literature for the shedding index-patient in this 
outbreak. However, researchers suggest that major 
outbreaks occurred because of a “super-spreader” 
index-patient having a disproportionately higher vi- 
ral load in their respiratory system (Xiao et al., 2018). 
Hence, the model was run applying a triangular dis- 
tribution based on the reported highest observed 
MERS-CoV concentration found in sputum of 
MERS-carrying patients (Min et al., 2016), using the 
reported minimum of 5.00 x 10°, average of 6.30 x 
10°, and maximum of 1.40 x 10’ copies/mL. Corman 
et al. (2015) even reported a concentration as high as 
10'' copies/mL, but this value was not used since re- 
searchers suggested it to be very unlikely, mention- 
ing that the probable highest concentration would 
more likely fall around 10’ copies/mL during the first 
week of shedding. Additionally, a two-day exposure 
duration (before diagnosis) was set based on the ob- 
served 44-hour exposure time in the South Korean 
MERS outbreak (Cho et al., 2016). Finally, a rate 
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of 3 ACH for the ventilation rate in the emergency 
room was used to match the outbreak conditions 
(Cho et al., 2016). Under these conditions, the cal- 
culated median risks of infection were found to be in 
the same ranges of the ones reported from Cho et al. 
(2016). Estimated median risks were found to be of 
2.08 x 10-7, 1.82 x 10-°, 9.82 x 107°, and 7.14 x 10° 
for nurses, healthcare workers, family visitors, and 
other patients sharing the room, respectively, while 
attack rates observed from the outbreak were mea- 
sured at 2.72 x 107%, 1.85 x 10-7, 5.56 x 10-7, and 
4.44 x 10-7, for same exposed groups. These compa- 
rable results provide validity to the model developed 
herein; thus, demonstrating its applicability to known 
outbreaks. 

Finally, the model developed estimated risk of 
infection from transmission to first-generation-only 
patients—those infected from direct exposure to in- 
dex patient—not from any potential second- or third- 
generation contacts. Therefore, risks might be under- 
estimated, as additional contacts would increase the 
exposure dose. 


5. CONCLUSIONS 


MERS, caused by the MERS-CoV virus, is be- 
lieved to have started in Arabian Peninsula. Due to 
the movement of people, the virus has a potential to 
cause MERS outbreak in other parts of the world, 
which was highlighted in a recent MERS outbreak in 
South Korea and China, and from the reported in- 
fected cases that occurred in 27 countries across all 
continents. In the main outbreaks, nearly all the cases 
were spread within a hospital setting, where the in- 
fected patient visited the hospitals, but was not quar- 
antined. In this context, this study used the QMRA 
approach to characterize the risk of MERS-CoV in- 
fection under hospital settings. The exposed popula- 
tions included in this study were the other patients 
sharing the same room, nurses, healthcare workers, 
and family visitors. 

The results showed that the nurses had the high- 
est daily risk of infection under a standard 6 ACH 
ventilation for typical hospital room, followed by 
healthcare workers, family visitors, while patients 
housed in the same room had the lowest daily risk 
of infection. Cumulative mean risk of infection was 
also highest for the nurses, which, by day 41, was 1.01, 
1.2, and 2.4 times the risk for the healthcare work- 
ers, family members, and other patients, respectively. 
Sensitivity analysis showed that the concentration of 
MERS-CoV in patient’s saliva was the most sensitive 
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parameter, constituting over 90% of the parameter 
uncertainty. The concentration of viruses in sputum 
is an inherently variable parameter that changes dur- 
ing the course of infection, illness, and recovery. Re- 
ducing the uncertainty contributed by this parameter 
is therefore unlikely. However, characterization of 
“super-spreader” in diagnostic laboratory tests could 
help prevent spread. 

Increasing the air exchange rate was found to be 
an effective risk reduction measure for the other pa- 
tients in the same room, but not for the other groups 
exposed to close-range airborne route. Using mask 
was found to be the most effective strategy, which 
could reduce over 90% of the risk for the exposed 
groups studied. Surgical and N95 masks has been re- 
ported to be highly effective (up to 100%) in prevent- 
ing transmission of respiratory diseases, however; 
some reports have found less than 50% of the health- 
care workers wearing masks even in the developed 
countries like the United States and Canada (Nichol 
et al., 2008; Park et al., 2004; Seto et al., 2003). 

The generalizable model developed herein us- 
ing the QMRA approach is intended to allow future 
risk assessors to adapt this framework to their spe- 
cific risk scenarios, by adapting each input parame- 
ter, accordingly based on newly available data. Such 
a model can be used to test hypotheses about control 
measures and risk management strategies. 
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